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Abstract With Teragauss magnetic fields, surface gravity sufficiently strong to signifi- 
cantly modify light paths, central densities higher than that of a standard nu- 
cleus, and rotation periods of only hundredths of a second, young neutron stars 
are sites of some of the most extreme physical conditions known in the Universe. 
They generate magnetic winds with particles that are accelerated to energies in 
excess of a TeV. These winds form synchrotron-emitting bubbles as the parti- 
cle stream is eventually decelerated to match the general expansion caused by 
the explosion that formed the neutron stars. The structure of these pulsar wind 
nebulae allow us to infer properties of the winds and the pulsating neutron stars 
themselves. The surfaces of the the stars radiate energy from the rapidly cooling 
interiors where the physical structure is basically unknown because of our im- 
precise knowledge of the strong interaction at ultrahigh densities. Here I present 
a summary of recent measurements that allow us to infer the birth properties of 
neutron stars and to probe the nature of their winds, the physics of their atmo- 
spheres, and the structure of their interiors. 
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Introduction 

Young neutron stars (NSs) probe some of the most extreme physical en- 
vironments in the Universe. Their rapid rotations and large magnetic fields 
combine to accelerate particles to extremely high energies, producing ener- 
getic winds that result in the slow spin-down of the stars and generate nebulae 
of synchrotron-emitting particles spiraling in a wound-up magnetic field. The 
structure of these nebulae is determined by the energy input from the central 
pulsars as well as the structure and content of the medium into which they ex- 
pand. In the centermost regions, relativistic outflows in the form of rings and 
jets are formed; the geometry of these emission regions reveals the orientation 
of the pulsar spin axes and can provide information on the formation of kicks 
imparted in the moments following their formation. Their large-scale struc- 
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tures reveal details of the magnetic field and signatures of interaction with the 
ejecta from the explosions that gave them birth. 

The stellar interiors are characterized by conditions and physical processes 
otherwise observed only within the nuclei of atoms. They are born hot, but 
cool rapidly due to neutrino production in their interiors. However, details 
of the interior structure of such stars remain poorly understood owing to our 
incomplete understanding of the strong interaction at ultrahigh densities and, 
since the neutrino production rate is critically dependent on the structure of the 
interior, the cooling rate is highly uncertain. In the standard cooling scenario, 
neutrino production proceeds primarily via the modified Urea process. Resid- 
ual heat diffuses from the core to the surface, manifesting itself as blackbody- 
like emission - modified by effects of any residual atmosphere - which peaks 
in the soft X-ray band. The rate at which the surface temperature declines 
depends critically upon the neutrino emission rate; thus, its measurement pro- 
vides constraints on hadronic physics at high densities. 

Perhaps the most stunning thing about NSs is the fact that we can actually 
make measurements that, directly or indirectly, probe the above properties. 
Particularly with the advent of sensitive high-resolution X-ray observations, 
we can now image jets and outflows from the wind termination shocks, iden- 
tify magnetic filaments in the nebular interiors, detect the thermal emission 
from shock-heated ejecta, and measure directly the pulsations from the rotat- 
ing stars and the emission from their ultra-hot surfaces. Here I describe the 
basic properties of young NSs and their nebulae, and summarize recent obser- 
vational work that has begun a revolution in our understanding of how these 
stars work. Brevity precludes a thorough review, and the reader is referred 
to recent articles by Kaspi et al., 2004 and Yakovlev and Pethick, 2004 for 
additional information and references. 

1. Pulsar Wind Nebulae 

Our basic understanding of PWNe stems from the picture presented by Rees 
and Gunn, 1974, and expanded upon by Kennel and Coroniti, 1984, in which 
an energetic axisymmetric wind is injected from a pulsar into its surroundings. 
As illustrated schematically in Figure 1 , the structure of a PWN is regulated 
by the input power from the pulsar and the density of the medium into which 
the nebula expands; the pulsar wind and wound-up toroidal magnetic field in- 
flates a bubble which is confined in the outer regions by the expanding shell of 
ejecta or interstellar material swept up by the SNR blast wave. The boundary 
condition established by the expansion at the nebula radius results in the 
formation of a wind termination shock at which the highly relativistic pulsar 
wind is decelerated to merge with the particle flow in the nebula. The shock 
forms at the radius r w at which the ram pressure of the wind is balanced by the 
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Figure 1. Schematic view of a pulsar and its wind nebula. See the text for a complete 
description. (Note the logarithmic size scaling in the PWN figure when comparing with images 
shown elsewhere in the text.) 

internal pressure of the PWN: 

r w = \/E/(4irr]cp), (1) 

where E is the rate at which the pulsar injects energy into the wind, 77 is the 
fraction of a spherical surface covered by the wind, and p is the total pres- 
sure outside the shock. Ultimately, the pressure in the nebula is believed to 
reach the equipartition value; a reasonable pressure estimate can be obtained 
by integrating the radio spectrum of the nebula, using standard synchrotron 
emission expressions, and assuming equipartition between particles and the 
magnetic field. Typical values yield termination shock radii of order 0. 1 pc, 
which yields an angular size of several arcsec at distances of a few kpc. 

As the relativistic fluid comprising the PWN encounters the freely-expanding 
ejecta, Rayleigh-Taylor instabilities result in the formation of a network of 
dense, optical line-emitting filaments (Jun, 1998). The density and magnetic 
field strength becomes enhanced in regions where the PWN encounters these 
filaments, producing enhanced synchrotron emission observed as radio fila- 
ments. Due to the pinching effect of the global toroidal magnetic field, the 
overall morphology of a young PWN is often elongated along the pulsar spin 
axis (Begelman & Li, 1992; van der Swaluw et al., 2004). Along the rota- 
tion axis the flow becomes collimated, producing jets. Pinch instabilities may 
disrupt the toroidal structure, however, changing the structure of the magnetic 
field in the outer nebula regions and relaxing the collimation of the jets far from 
the pulsar (Begelman, 1998). 

The overall geometry of the PWN, as well as that of the emission from jets 
or ring-like structures near the termination shock, thus provides a direct indica- 
tion of the pulsar geometry. The details of the jet morphology and the emission 
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Figure 2. Chandra images of the Crab Nebula (left), G54. 1+0.3 (center), and PSR B 1509-58 
(right) showing the complex emission from these PWNe, including jet outflows and toroidal 
structures. 



structure in the postshock region provide the strongest constraints available on 
wind composition and particle acceleration in PWNe. For cases in which the 
pulsar proper motion is also known, constraints on the kick velocity mecha- 
nism can be derived based on the degree of alignment between the velocity 
vector and the pulsar spin axis. In later stages the PWN interacts with the 
reverse shock formed in the SNR in which the NS was born. This interac- 
tion causes the disruption of the PWN, often leading to composite SNRs with 
complicated PWN structures in their interiors. 

1.1 Jets and Tori 

In the inner portions of the Crab Nebula, optical wisps mark the position 
of the wind termination shock, at a distance of ~ 0.1 pc from the pulsar. The 
brightness and position of these wisps varies in time, with inferred outflow 
speeds up to 0.7c (Hester, 1998). As shown in Figure 2 (left), high resolution 
X-ray images reveal a ring of emission at the position of the wisps (Weisskopf 
et al, 2000), providing a direct connection between the unshocked pulsar wind 
and the bulk properties of the nebula. Material from the inner ring forms a 
series of toroidal X-ray wisps that are variable with time (Hester et al., 2002). 
The geometry of these X-ray features imply a tilted torus, and a jet of ma- 
terial flows perpendicular to the plane of the toroid, extending some 0.25 pc 
from the pulsar. A faint counterjet is also observed, along with significantly 
enhanced X-ray emission from the leading portion of the toroid, presumably 
the result of Doppler beaming. One troubling aspect of this suggestion is that 
the brightness distribution around the inner ring does not match that of the 
outer toroid; indeed, the brightness is rather uniform except for some small 
clump-like structures that vary in position and brightness with time. 

A handful of other PWNe display X-ray features that suggest the presence 
of extended ring-like structures and narrow collimated components. The size 
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of the ring-like features places the emission region near the pulsar wind ter- 
mination shock. The spectral and temporal properties of the collimated struc- 
tures argue that they are focused jets of high speed material, as observed in 
the Crab. Such observations have already begun to inspire new axisymmet- 
ric MHD models that predict similar features (e.g., Komissarov & Lyubarsky 
2004), and ongoing observational studies promise to further constrain and re- 
fine such models. In particular, the confining mechanism for jets is not well- 
understood; many jets display some amount of curvature, with the Vela pulsar 
jet being an extreme example in which the morphology is observed to change 
on timescales of months (Pavlov et al, 2003). This may be the result of pinch 
instabilities disrupting the toroidal structure of the confining magnetic field 
(Begelman 1998), or could be indicative of an interaction of the jet material 
with the ambient medium. There also appears to be a wide variation in the 
fraction of spin-down energy channeled into the jets, ranging from roughly 
2.5 x 10" 5 for PSR J0205+6449 in 3C 58 to nearly 10~ 3 for PSR B 1509-58. 
And, while Doppler beaming is invoked to explain the large brightness vari- 
ations in jets and the associated counterjets, as well as around the observed 
toroidal structures, it is not clear that this alone is sufficient to explain the ob- 
servations. 

Chandra observations of G54. 1+0.3 (Lu et al, 2002) reveal a central 136 ms 
pulsar (Camilo et al., 2002a) embedded in a diffuse lf5 x 1'.2 nebula (Figure 
2, center). The pulsar is surrounded by an X-ray ring for which the X-ray 
emission is brightest along the eastern limb. When interpreted as the result of 
Doppler boosting, this implies a post-shock velocity of ~ 0.6c (Lu et al. 2002; 
Romani & Ng 2003). Faint bipolar elongations running roughly east-west, 
perpendicular to the long axis of the ring, are also observed. These apparent 
outflows, which presumably lie along the pulsar rotation axis, are more diffuse 
than the jets in the Crab Nebula, yet appear to carry away a considerably larger 
fraction of the energy; they comprise roughly the same luminosity as the cen- 
tral ring, which is in stark contrast to the Crab where the torus outshines the 
jets by a large factor. 

Chandra observations of PSR B1509— 58 (Gaensler et al, 2002) demon- 
strate that this young and energetic pulsar associated with G320.4-1.2 powers 
an extended and extremely complicated PWN, with structures on scales from 
~ 10' down to the spatial resolution limit (Figure 2, right). The elongated 
PWN has a clear axis of symmetry centered on the pulsar, presumably repre- 
senting the projected orientation of the pulsar spin axis. To the southeast of 
the pulsar, the nebula is dominated by a narrow jet-like feature approximately 
6 pc in length. The lack of a similar feature to the north can be explained 
by Doppler boosting if the pulsar's spin axis is inclined to the line-of-sight by 
£ 30° (Gaensler et al., 2002). In the central regions of the PWN, a pair of semi- 
circular arcs lie ~ 0.5 and ~ 1 pc to the north of the pulsar. Gaensler et al. 
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Figure 3. Left: Chandra image of 3C 58. Complex filamentary loops fill the interior region. 
Right: The innermost region of 3C 58 showing the NS embedded in an elongated structure. A 
curved jet extends to the west, with a hint of a counterjet component in the east. 



(2002) note that if the inner region of these arcs represents the position of the 
pulsar wind termination shock, then the flow time to the arcs is much shorter 
than the synchrotron lifetime of the emitting particles based on equipartition 
estimates of the magnetic field. Thus, unlike for the Crab torus, where these 
timescales are similar, the emission from the arcs is not the result of large syn- 
chrotron cooling at this position. Instead, the arcs appear to resemble the series 
of concentric wisps seen for the Crab which are interpreted as sites of electron 
compression in an ion-dominated flow (Gallant and Arons, 1994, Gaensler et 
al., 2002). 

The innermost region of 3C 58 (see Figure 3) consists of a bright, elongated 
compact structure centered on the pulsar J0205+6449. This inner nebulosity 
is bounded along the western edge by a radio wisp (Frail and Moffett, 1993), 
and is suggestive of a toroidal structure that is tilted about a north-south axis, 
with the pulsar at its center. The eastern side of the toroid is slightly brighter 
than the western side, suggesting that the eastern side is beamed toward us. If 
interpreted as a circular termination shock zone, the inferred inclination angle 
in the plane of the sky is roughly 70 degrees (Slane et al., 2002). 

The elongated structure extending westward from the position of the pulsar 
has the appearance of a jet (Figure 3, right). Its orientation is consistent, in 
projection, with the pulsar rotation axis inferred from the wind termination 
shock region discussed above, and also the east-west elongation of the entire 
PWN (Figure 3, left). The structure shows considerable curvature, similar to 
that seen in the Crab Pulsar jet. A faint structure that may be a counterjet is 
observed to the east of the pulsar. The observed luminosity is nearly a factor 
of 10 smaller than that for the torus. For the Crab Nebula, the torus is nearly 
20 times more luminous than the jet in X-rays, while for PSR B 1509-58 the 
jet is brighter than the extended inner emission (Gaensler et al. 2002). 

The jet/torus morphology observed in these PWNe provides the geometry of 
the pulsar system, yielding both the projected direction of the spin axis and the 
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inclination angle. Modeling of such emission in other PWNe holds promise for 
understanding the kicks that give pulsars their large space velocities (Ng and 
Romani, 2004). The jets observed in the Crab and Vela pulsars, for example, 
are aligned with their proper motion vectors (Aschenbach & Brinkman 1975; 
Helfand et al. 2001). If the kick that gave these pulsars their proper motion 
was generated in the supernova explosion by some asymmetric mass ejection, 
then this alignment requires an initial pulsar spin period that is short relative 
to the kick timescale, so that the impulse of the kick is averaged over many 
rotations of the star (Lai et al., 2001). Romani and Ng, 2003 reach similar 
conclusions for PSR J0538+2817 in the supernova remnant (SNR) S147. By 
modeling the faint extended PWN emission as a jet and torus, they derive a 
spin axis direction that is aligned with the vector from the SNR center to the 
current pulsar position. For some pulsars [e.g. J0205+6449 in 3C 58 (Murray 
et al., 2002) and J1811-1925 in Gll.2-0.3 (Kaspi et al., 2001)], we believe 
that the initial spin period was much longer than typical pulsar kick timescales. 
This would suggest that their proper motions should not necessarily be aligned 
with the jet direction. Future radio timing observations of these pulsars will 
ultimately lead to such proper motion measurements. 

1.2 Filaments in PWNe 

Extensive filamentary structure is observed in Ha, [OIII], and other optical 
line images of the Crab Nebula. Based on their observed velocities, these fila- 
ments form an expanding shell of ejecta that surrounds the nonthermal optical 
emission from the nebula. High resolution images with HST reveal detailed 
morphology and ionization structure suggesting that the filaments form from 
Rayleigh-Taylor instabilities as the expanding relativistic bubble encounters 
slower moving ejecta (Hester et al., 1996), a picture supported by MHD sim- 
ulations that show that 60-75% of the swept-up mass ends up concentrated in 
such filaments (Jun 1998, Bucciantini et al. 2004). Radio observations reveal 
filaments that coincide with these optical filaments, presumably corresponding 
to synchrotron emission from regions of enhanced density and magnetic field 
in the form of magnetic sheaths that form as the pulsar-injected energy en- 
counters the thermal filaments (Reynolds, 1988). Such filamentary structure is 
not observed in X-rays, however, suggesting that the electrons with sufficient 
energy to radiate X-rays do not reach the shell of filaments. This is consistent 
with the observed smaller extent of the X-ray emission in the Crab nebula rel- 
ative to its radio size, and indicates a larger magnetic field than is observed in 
3C 58 and PSR B 1509-58. 

Recent Chandra observations of 3C 58 reveal a complex of loop-like fila- 
ments most prominent near the central regions of the PWN (Figure 3, left), but 
evident throughout the nebula (Slane et al., 2004a). These structures, whose 
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X-ray spectra are nonthermal, are very well correlated with features observed 
in the radio band (Reynolds and Aller, 1988). Optical observations reveal faint 
thermal filaments as well (van den Bergh, 1978), which presumably have an 
origin similar to that of the Crab filaments. The velocities of these optical fila- 
ments in 3C 58 are ~ ±900 km s _1 (Fesen, 1983), sufficiently high to indicate 
that the PWN is young, but too small to account for the current size of 3C 58 
if the historical age is assumed - one of several standing problems with regard 
to its evolution (Chevalier, 2004). A detailed comparison of the X-ray and op- 
tical images shows that most of the X-ray filaments do not have corresponding 
optical structures, however. While comparisons with deeper optical images are 
clearly needed, the fact that many of the X-ray features without optical coun- 
terparts are brighter than average in X-rays suggests that these may actually 
arise from a different mechanism. Slane et al. (2004) propose that the bulk of 
the discrete structures seen in the X-ray and radio images of 3C 58 are mag- 
netic loops torn from the toroidal field by kink instabilities. In the inner nebula, 
the loop sizes are similar to the size of the termination shock radius, as sug- 
gested by Begelman (1998). As the structures expand, they enlarge slightly as 
a consequence of the decreasing pressure in the nebula. Some of the observed 
X-ray structure in the outermost regions may be the result of thermal filaments 
produced by Rayleigh-Taylor instabilities, similar to the filaments in the Crab 
Nebula. A shell of thermal X-ray emission demonstrates the presence of ejecta 
in these outer regions (Bocchino et al., 2001, Slane et al., 2004a. 

It is worth noting that considerable loop-like filamentary structure is evident 
in Chandra observations of the Crab Nebula as well (Weisskopf et al., 2000). 
These features are primarily observed encircling the bright Crab torus, per- 
pendicular to the toroidal plane, and may result from currents within the torus 
itself. It is at least conceivable that such currents are signatures of the kink 
instabilities suggested above. 

1.3 Large-Scale Structure of PWNe 

The large-scale elongated shape of 3C 58 is similar to that found (particu- 
larly in the radio band) for a number of other PWNe including the Crab Neb- 
ula and G54. 1+0.3. Magnetohydrodynamical calculations by Begelman & Li, 
1992, and van der Swaluw, 2004, show that such an elongation can result from 
the pinching effect of a toroidal magnetic field for which the projected axis lies 
along the long axis of the PWN. The pinching effect results in a low pressure at 
the edge of the bubble along the major axis with respect to the (higher) pressure 
at the edge of the minor axis, which yields the elongated structure. The elon- 
gation thus marks the projection of the spin axis of the pulsar producing the 
wound-up field. In 3C 58 this is consistent with the inference of an east-west 
direction for the projected spin axis based on the interpretation of the extended 
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Figure 4. Left: Chandra image of G292.0+1.8. The inset shows the central region at energies 
above 4 keV, where the pulsar and its wind nebula dominate. Right: Schematic diagram of the 
evolutionary state of G292.0+1.8. The lack of iron observed in the spectrum indicates that the 
reverse shock has not yet made its way to the center of the remnant, where the PWN resides. 



structure in the inner nebula as being associated with a tilted ring-shaped wind 
termination shock zone (Slane et al., 2002). 

The structure of a PWN can be altered significantly through interaction with 
the reverse shock from the SNR in which it resides. In its early evolution the 
PWN is basically freely-expanding, encountering only small amounts of slow- 
moving ejecta in the SNR interior. As the SNR blast wave sweeps up sufficient 
amounts of circumstellar/interstellar material, a reverse shock is driven back 
through the ejecta. As this reverse shock propagates, heating the ejecta, it 
will eventually reach the PWN. Chandra studies of the oxygen-rich remnant 
G292.0+1.8 reveal an SNR in the intermediate stages of this process. The 0.5- 
10 keV X-ray image is presented in Figure 4 (left), and shows the complex 
structure associated with the shock-heated ejecta and CSM (Park et al., 2002). 
The inset shows the central image at energies above 4 keV, and reveals a com- 
pact pulsar surrounded by a wind nebula (Hughes et al., 2001, Camilo et al., 
2002b, Hughes et al, 2003). X-ray spectra of the SNR show metal-rich ejecta 
with strong lines of oxygen and neon, but a distinct shortage of iron emission 
(Park et al., 2004), indicating that the reverse shock has not yet propagated 
sufficiently far toward the center to heat the iron-rich material that was formed 
closest to the core of the progenitor (Figure 4, right). 

The morphology of the PWNe in G327. 1 — 1.1 presents a rather different 
picture. The radio image (Figure 5, left) reveals a well-defined SNR shell with 
a bright PWN in its interior, distinctly offset from the geometric center. A 
finger of emission extends to the northwest of the radio PWN, and Chandra 
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Figure 5. Left: MOST image of the composite SNR G327. 1 — 1.1. Right: Chandra image 
of the PWN in G327. 1 — 1.1 (with radio contours), showing diffuse emission surrounding a 
compact, but resolved, X-ray source (inset). 



observations reveal a compact X-ray source at this location (Figure 5, right). 
The source is slightly extended (see inset) suggesting that we are seeing ma- 
terial near the wind termination shock. The compact source resides at the tip 
of a trail of emission that leads back to the bulk of the radio nebula, rather 
than at the center of the PWN (as in G292.0+1.8, for example), suggesting that 
the PWN morphology results from a combination of the pulsar motion and the 
passage of the reverse shock which has apparently disrupted the western side 
of the nebula (Slane et al., in preparation). 

2. Neutron Star Cooling 

The cooling rate of isolated NSs has been a subject of considerable theoreti- 
cal work predating even the discovery of the first pulsars (e.g., Bahcall & Wolf 
1965). The poorly understood properties of the strong nuclear potential at the 
densities found in NS interiors make these calculations difficult, and lead to a 
wide range of predictions based on different assumptions for the equation of 
state, composition, and details of superconductivity (see, e.g., reviews by Tsu- 
ruta, 1998, and Yakovlev & Pethick, 2004). While there is a clear consensus 
that the early cooling proceeds via neutrino emission from the NS core, the 
timescale over which this dominates depends critically on the neutrino produc- 
tion rate which, in turn, can vary by orders of magnitude depending upon the 
state of matter in the interior. 

Broadly speaking, models can be divided into "standard" and "non-standard" 
cooling scenarios based on the rate of neutrino production in the NS interior. 
At the high densities and low proton fractions expected in NS interiors, the 
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Figure 6. Left: Representative neutron star cooling curves for neutrino rates corresponding to 
different interior conditions. Data from Page, 1998. Right: Temperature vs. Age measurements 
for neutron stars, along with plots for standard cooling (bold), "minimal" cooling, and mass- 
dependent cooling invoking dUrca (see text). 

direct Urea reactions (n — > p + e _ + v e and p + e~ — > n + v e ) cannot conserve 
both energy and momentum. Instead, a bystander baryon is required for each 
interaction to absorb momentum. The neutrino rate for this so-called modified 
Urea (mUrca) process is considerably lower than in the direct (dUrca) process 
because of the extra interaction required, and is the basis for standard cooling 
models, neutrino bremsstrahlung, and plasmon neutrino processes). In Figure 
6 (left), we plot cooling curves for different models of the NS interior and its 
properties using representative neutrino rates (see Page, 1998 and references 
therein). The solid curve corresponds to "standard" cooling using an equation 
of state of moderate stiffness. 

Non-standard cooling models incorporate neutrino emissivities associated 
with other processes that may operate in NS interiors, such as those arising 
from the presence of pion condensates which may form at sufficiently high 
densities. The resulting pion-induced beta decay leads to very a high neutrino 
emissivity and a correspondingly shorter cooling time for the NS interior. Sim- 
ilar processes involving kaon condensates or quark matter may operate as well. 
Alternatively, equations of state that allow a high proton fraction in the interior 
may allow the dUrca process to proceed; this also leads to extremely high neu- 
trino production rates (e.g., Kaminker et al., 2002). These nonstandard cooling 
mechanisms modify the NS cooling curves substantially. The dashed curves in 
Figure 6 represent approximations for several nonstandard cooling models and 
illustrate the associated rapid cooling (Page 1998). 

The effects of superfluidity can substantially moderate such rapid cooling 
because the significantly reduced heat capacity of the superfluid particles re- 
duces the neutrino rate considerably. Thus, slow cooling can occur from a com- 




12 



bination of slow neutrino rates and a high degree of superfluidity, while rapid 
cooling places strong constraints on both, and thus on the structure and physics 
of the NS interior, as well as on the microphysics of superfluidity. In partic- 
ular, as the temperature approaches the critical temperature for the superfluid 
state, the formation and breaking of Cooper pairs opens another channel for 
neutrino emission that can lead to more rapid cooling. Page et al., 2004 have 
calculated a "minimal cooling" scenario which extends the standard mUrca 
cooling scenario to include the contributions from this Cooper pair process. 
As discussed below, for some superfluidity models they find sufficiently rapid 
cooling to explain most observations of young NSs. Yakovlev and Pethick, 
2004, on the other hand, have considered both standard and enhanced cooling 
models. Using the NS mass as a free parameter, they find models for which 
the dUrca process becomes active for sufficiently high masses, thus leading to 
a picture in which the young, cool NSs correspond to those with higher mass 
(see below). 

2.1 Measuring Neutron Star Temperatures 

Treating the emission from the NS surface as a blackbody, X-ray spectral 
fitting provides a measure of the gravitationally redshifted temperature and 
luminosity (assuming the distance is known): 



rpOO 
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where the quantities on the left are observed at infinity, and those on the right 
are at the NS surface; M and R are the mass and radius of the NS. From 
this we calculate the effective radius, which can be compared directly with 
predications for different equations of state: ^// = L°° /[4ira(T^) 4 } (where 
a is the Stefan-Boltzmann constant). Alternatively, if only an upper limit on 
the source luminosity is determined, a temperature upper limit can be derived 
by assuming a value for the NS radius. 

As with all stars, the emission from the surface of a NS is not a blackbody; 
rather, it is modified by the presence of whatever atmosphere might exist. One 
expects the surface of the NS to be covered with Fe, but an atmosphere con- 
sisting of H, He, and/or intermediate-mass elements acquired either from ejecta 
fallback following the neutron star's formation, or from material accreted from 
the ISM, is also a possibility. From models of nonmagnetic atmospheres, the 
primary effect of H or He atmospheres is a considerable deviation of the high 
energy end of the spectrum relative to the Wien tail of a pure blackbody. The 
result is that attempts to fit the observed emission with a blackbody model will 
overestimate the effective temperature - typically by as much as a factor of 
two. The inferred size of the NS would, in turn, be underestimated in order to 
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yield the same flux. For atmospheres dominated by heavier elements the effect 
is considerably reduced, and the blackbody fit gives a good approximation to 
the temperature. 

2.2 Confronting Cooling Models 

In Figure 6 (right) we plot the measured temperatures, or upper limits, for 
the pulsars and compact objects in SNRs for which these values are best deter- 
mined. Values for known pulsars are plotted with closed circles. For compar- 
ison, curves are shown for standard cooling as well as the "minimal" cooling 
model of Page et al., 2004 (dashed curve). The latter model assumes no en- 
hanced cooling mechanisms in the interior. The primary mechanism which 
results in faster cooling is neutrino emission from the Cooper pair breaking 
and formation process, which is heavily dependent upon the assumed super- 
fluidity model (here we have plotted their model "a" for the neutron 3 P2 gap); 
the predicted cooling rate is sufficiently rapid to explain the measurements for 
most pulsars, although the temperatures for the Vela Pulsar and J0205+6449 
in 3C 58 fall appear to require more rapid cooling. The inferred temperature 
for RX J0007. 0+7302 (plotted as an open box in Figure 6), a compact X-ray 
source in the SNR CTA 1, also falls below this minimal cooling scenario (Slane 
et al., 2004). While pulsations have yet to be detected from this source, Chan- 
dra observations by Halpern et al., 2004 reveal a slightly extended source with 
a distinct jet-like tail that is the signature of the pulsar outflows described in 
Section 1.2. 

For comparison, cooling curves illustrating the effects of increasing mass 
are plotted as solid curves, using the model lp for proton superfluidity from 
Yakovlev and Pethick, 2004. Also plotted in the Figure, as dashed arrows, 
are upper limits based on a neutron star search in nearby SNRs (Kaplan et al., 
2004). In this study, deep X-ray observations were used to identify, sources 
within the field of the SNRs, and follow-up optical and IR observations were 
used to screen for non-NS counterparts. While the lack of detection of a NS 
in these remnants (G127. 1+0.5, G84.2-0.8, G93.3+6.9, and G3 15.4-2.3 in in- 
creasing order of temperature upper limit) may indicate that these SNRs origi- 
nated in Type la events, or all formed black holes, this is statistically unlikely. 
Future measurements from this study will solidify this picture, and perhaps 
provide further evidence for rapid cooling in young neutron stars. 

3. Summary 

Due in large part to the availability of high resolution X-ray measurements, 
the study of young NSs and their PWNe have yielded dramatic new informa- 
tion on the nature of the stellar interiors and the structure of their winds. It 
is now clear that the axisymmetric wind from a pulsar goes through a termi- 
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nation shock as it joins the slower flow of its extended nebula, and that jets 
and toroidal structures characterizing the inner structure can be used to infer 
the orientation of the pulsar spin axis. The brightness variations in the inner 
nebula and jets, as well as the spectral and spatial structure of these regions, 
yield information on the outflow geometry as well as the fraction of spin-down 
energy being channeled into these regions. With broad application to particle 
acceleration and jet formation in astrophysical settings, these observations are 
providing constraints on theoretical models of considerable importance for a 
wide range of problems. Similarly, new observations are providing unprece- 
dented capabilities for detecting young neutron stars in SNRs, and for char- 
acterizing their emission. These have forced a revised look at models for the 
structure of NS interiors and provide the best opportunity for addressing the 
possibility that exotic states of matter reside in NS cores. Through additional 
and more sensitive observations of these systems, we anticipate significant re- 
finements, and undoubtedly new surprises, in broad picture of their structure 
that is currently unfolding. 
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